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ABSTRACT 


Context. Co-orbital systems are bodies that share the same mean orbit. They can be divided into different families according to the 
relative mass of the co-orbital partners and the particularities of their movement. Janus and Epimetheus are unique in that they are the 
only known co-orbital pair of comparable masses and thus the only known system in mutual horseshoe orbit. 

Aims. We aim to establish whether the Janus-Epimetheus system might have formed by disruption of an object in the current orbit of 
Epimetheus. 

Methods. We assumed that four large main fragments were formed and neglected smaller fragments. We used numerical integration 
of the full N-body problem to study the evolution of different fragment arrangements. Collisions were assumed to result in perfectly 
inelastic merging of bodies. We statistically analysed the outcome of these simulations to infer whether co-orbital systems might have 
formed from the chosen initial conditions. 

Results. Depending on the range of initial conditions, up to 9% of the simulations evolve into co-orbital systems. Initial velocities 
around the escape velocity of Janus yield the highest formation probability. Analysis of the evolution shows that all co-orbital systems 
are produced via secondary collisions. The velocity of these collisions needs to be low enough that the fragments can merge and not 
be destroyed. Generally, collisions are found to be faster than an approximate cut-off velocity threshold. However, given a sufficiently 
low initial velocity, up to 15% of collisions is expected to result in merging. Hence, the results of this study show that the considered 
formation scenario is viable. 


1. Introduction 


In a co-orbital system, two or more bodies share a mean orbit 
around a central body in a 1:1 orbital resonance. The possible or¬ 
bits have been described by Lagrange in his Analytical Mechan¬ 
ics in 1788. We speak of tadpole orbits when one body librates 
around the Lagrange points L4 or L5 of another. Libration around 
L3, L4, and L5 constitutes a horseshoe orbit ( |Brown|191 l| l. Many 
trojan asteroids are in co-orbital trajectories with planets, but 
only in the Saturnian system are satellites known to be in co¬ 
orbital motion. While there are many examples of bodies with 
much smaller co-orbital partners, the Saturn moons Janus and 
Epimetheus are a unique system of two bodies with similar mass 
in horseshoe orbit. This co-orbital motion was confirmed in 1980 
by the Voyager I probe (Synnott et al.|198T]l. Harrington & Sei 


delmann| ( |1981| l showed the stability of this system through nu 
merical integrations including the oblateness of Saturn. [Dermott| 
& Murra^ l 1981[ ) developed an analytical study estimating the 
radial and angular amplitudes of each horseshoe orbit. A theory 
of motion was proposed two years later by | Yoder et’ak] ( |1983| l, 
who explored the phase space of co-orbital systems with compa¬ 
rable masses, identifying the equilibrium points and showing the 
possible trajectory types (tadpole and horseshoe). Considering a 


system with up to nine equal-mass co-orbital bodies, Salo & Yo 
|der|( 1988]l found the location of the equilibrium points and their 


stability. Eollowing a Hamiltonian approach, |Sicardy & Dubois] 
( 2003| l demonstrated the connection between co-orbital systems 
of comparable masses with that of the restricted case. 

The formation of co-orbital systems has been the subject of 
intensive studies. In 1979, Yoder ( |1979 1 proposed three possible 
mechanisms for the formation of the Jovian trojans: 


1. early capture during the accretion of Jupiter, 

2. the continuous capture of asteroids, and 


3. transfer of Jovian satellites to the Lagrange points. 


In the first scenario, gas drag plays an important role. Peale 
( |1993| l and |Murray| ( |1994) l showed that, in the presence of a 
gas, L4 and L5 can still be stable for proto-Jupiter. In 1995, 
Kary & Lissauer|( 19951 showed that a protoplanet on an eccen¬ 


tric orbit can indeed trap planetesimals in a 1:1 resonance, de¬ 
caying from a large-libration horseshoe orbit to a tadpole orbit 
around L5. Eleming & Hamilton ( 2000) l explored the effect due 
to mass accretion and orbital migration on the dynamics of co¬ 
orbital systems. They showed how these perturbations affect the 
radial and angular amplitudes of the trajectories. [Kortenkamp 
( |2005| l later proposed a mechanism to efficiently capture plan¬ 
etesimals on a satellite orbit after gravitational scattering by the 
protoplanet, making this a more likely fate than trapping in a 
co-orbital. IChanut et al.| ( |2008l l determined the size limit of trap¬ 
ping planetesimals via nebular gas drag to be 0.2m per AU of 
orbit of the protoplanet. They ( Chanut et al.|2013'] l later extended 
their approach to a non-uniform gas with a cavity in the vicin¬ 
ity of a heavy planet, but found little to no co-orbitals for ec¬ 
centricities comparable to Janus and Epimetheus. All these sce¬ 
narios result in co-orbital configurations with a large secondary 
body and a small partner, that is, with a low mass ratio. This al¬ 
most excludes trapping as a mechanism for the formation of the 
Janus-Epimetheus system. Eurthermore, |Hadjidemetriou & Voy^ 
|atzis| ( |20TT| l showed that a system of two equal-mass bodies in a 
1:1 resonance tends to evolve under gas drag into a closely bound 
binary satellite system. Another take on the subject by |Cresswell| 
|& Nelson] ( |2006| ) was the study of the evolution of a swarm of 
similarly sized planetesimals in a gas disk. They showed that the 
formation of mean-motion resonances, including 1:1 resonances, 
is commonplace, but that the bodies are still prone to be lost to 
accretion by the central body. Giuppone et al. P012| l followed 
a similar approach to [Chanut et al.| ( |2013| l, but studying a hypo- 
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thetical extrasolar system. They showed that a pair of inwardly 
migrating planets can enter into a 1 ; 1 resonance at the edge of a 
cavity in a gas disk. Whether their result is applicable to Janus 
and Epimetheus is not known to us. 

Laughlin & Chambers (|2002 1 showed that stable 1:1 reso¬ 


nances of equal-sized planets around a central body exist. They 
hypothesized that at the L 5 point of a planet in the accretion disk 
of a protostar, another planet could accrete. Beauge et al. p007| l 
investigated formation by accretion of a swarm of planetesimals 
in a tadpole orbit around L 4 or L 5 of a giant extrasolar planet. 
They showed that an Earth-type body can form, but with a very 
low efficiency and thus a low final mass. Chiang & Lithwick 
P005| l considered this a likely process for the formation of the 
Neptune trojans. Izidoro et al. p010| l studied the mechanism for 
the Saturnian system and showed that formation of co-orbital 
satellite pairs with a relative mass of 10 “'^ to 10 “® is possible. 
This again rules out the formation of a co-orbital pair with simi¬ 
lar masses. 

While much research has been done on the subject of co¬ 
orbital systems, there is as yet no satisfactory explanation for 
the origin of the Janus-Epimetheus system in particular. Our 
approach here is based on the collision of two bodies at or 
close to the current orbit of Epimetheus. This is consistent with 
the equal densities (see Table [T]) and compositions of Janus 
and Epimetheus. Moreover, results by |Charnoz et al.| p010| l 
show a recent formation of the Saturn moonlets from its rings 
through series of subsequent collisions. [Crida & Charnoz]p012| l 
later extended this scenario to all regular moons of Saturn. In 
this context, low-velocity collisions are frequent and are al¬ 
most inevitable, especially for the predecessors of Janus and 
Epimemetheus. We assume the disruption into four main frag¬ 
ments that later collide again constructively to form a pair of 
bodies in a co-orbital configuration. We establish the statistical 
likelihood of this outcome in 8000 numerical simulations. We 


used Chambers’ Mercury (Chambers 19991 N-Body integrator 
with a Burlisch-Stoer algorithm. 

In Sect. 2, we explain our model in detail and introduce the 
methods used to obtain our results. Section 3 shows the setup 
and results of our simulations. Einally, we discuss the viability 
of our model considering cases with six and eight fragments, and 
propose opportunities for future studies in Sect. 4. 


2. Model and method 

We investigated whether the Janus/Epimetheus co-orbital system 
might have formed by disruption of a large body in the mean or¬ 
bit of Epimetheus. Eigure[T][a,b] illustrates the scenario that we 
believe can lead to co-orbital configurations. We assumed that 
four large fragments are formed by a collision event. While the 
initial collision is not subject of this study, we note that a low col¬ 
lision velocity and low density (compare Table [T} of the collid¬ 
ing bodies is a requirement. The collision is neither completely 
destructive nor is it constructive. After the collision event, frag¬ 
ments spread out on orbits around Saturn. Fragments can be lost 
by ejection from the system, collision with Saturn, or collision 
with each other. In the last case, fragments can unite to form a 


larger body, given a sufficiently low collision velocity (see Ko- 
[rtenkamp & Wetherill|2000| ). 

Our simulations start right after the initial collision event, 
with the four fragments a, ff, y, and 6 on the vertices of a square, 
moving in opposite directions with equal speed (see Fig. I^for 
details). We always chose equal masses for bodies opposite on 
the square. The length ao of the sides of the square was chosen 
to be flo = + do , where Ri is the radius of the largest body 
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Fig. 2: Initial configuration of the simulations, v, are given in 
the reference frame moving with vq. The initial velocity of each 
fragment in the simulation frame of reference is given by v,- = 
V, -H vq. The direction of v, is given by the vector from the centre 
of the square to the respective vertex. The collision fragments a, 
P, y, and 6 are positioned on the vertices of a square centred on 
To, the position vector with respect to Saturn. The angle ip gives 
the rotation of the square around the vector Tq x vq (the normal 
of the orbital plane). For cp — 0 a, y, and Saturn are connected 
by a straight line, with y between a and Saturn. The fragments 
are always set on the vertices in a counter-clockwise sense, as 
shown in the sketch. The dimensions here are not to scale. 


and do - 100m. The parameters varied in this study include the 
masses m, of the fragments (given as fractions of the mass of 
Janus and Epimetheus, see Tables [T] and |^, the orientation of 
the square in the orbital plane denoted by the angle <p (going 
from 0 to TT because of symmetry), and the initial speed Uinit (see 
Fig. 1 ^ of the fragments in the coordinate system moving with 
Vq relative to Saturn. The position and velocity {ro,Vo) of the 
centre of the square are equal to the state vector of Epimetheus 
at an arbitrarily chosen point of its orbit. This is the same for all 
simulations. 

We used the package Mercury to numerically integrate the 
system using a Burlisch-Stoer algorithm ( [Chambers 1999[ ). We 
included the moons Mimas and Tethys into our simulations, 
since they exert the strongest influence on the Janus-Epimetheus 
system, but neglected other objects. We also considered the 
oblateness of Saturn through J 2 , J 4 , and Je (Table |^. To obtain 
the state vectors of Mimas, Tethys, and Epimetheus (centre of 
the square from which the fragments start), we used the data 
given by [Jacobson et al. ( 2008|l and followed the approach pre¬ 
sented by Renner & Sicardy ( 2006[ l to transform the geometric 
elements. We took the mean distance between the system Janus- 
Epimetheus and Saturn as the unit of length Rq. The simulation 
code was set to record the state vectors of all bodies every ’/loo 
day. We need this high time resolution to study the velocities 
and angles of collisional events between fragments. Mercury as¬ 
sumes a perfectly inelastic collision where the colliding bodies 
aggregate into a new body of the added masses of the colliding 
bodies, and the momentum is preserved. We therefore needed to 
analyse the collisions to ensure that this assumption is realistic. 
The outcome of a collision is strongly dependent on the colli¬ 
sion velocity and the physical characteristics (density, strength, 
etc.) of each body. [Kortenkamp & Wetherill| ( [2000] l established 
an escape velocity term 


2G (mi + m 2 ) 
Ri + R2 


( 1 ) 
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(a) Collision event. Body A collides with body B on (b) Disruption into four main fragments. The remaining 
Epimetheus’ orbit. debris is assumed to be negligible. 




(c) The main fragments distributed on orbit. 


(d) Eragments collide to form Epimetheus. 



(e) Fragments collide to form Janus (f) Janus and Epimetheus enter into horseshoe orbits. 

Fig. 1: Sketch of typical co-orbital formation scenario from our simulations. 


for the encounter of rocky bodies, where m,- is the mass of body i, 
Ri is its radius, and G is the gravitational constant. For collision 
velocities below Umax, a collision can be assumed to be construc¬ 
tive. For a collision of Janus and Epimetheus, Umax ~ 58m/s. 
However, their study involved bodies with a density about five 
times as high as that of Janus and Epimetheus. Since in a colli¬ 
sion event between fluffy bodies (such as Janus and Epimetheus, 
whose density is lower than 700kg/m^) kinetic energy is lost to 
deformation, higher velocity collision events can still be con¬ 
structive. High-velocity collisions can also result in the catas¬ 
trophic disruption of one or both colliding bodies. Multiple stud¬ 
ies on the disruption threshold (|Colwell|1994 Benz & Asphaug 


1999 1 Leinhardt & Stewart [2009 |l have been conducted. Using 

Fig. 11 from Leinhardt & Stewm^ ( |2009| ) and a target radius of 
6 ■ 10"^m (the rough size of Epimetheus), we obtain a disruption 


threshold of Q'^ =» 6.5 ■ lO^J/kg. This gives a cutoff velocity 
UD = J2Q0— (2) 

Y ifip 

(where nij is the mass of the target and nip the mass of the projec¬ 
tile) of approximately 60m/s for a collision with Janus. There¬ 
fore we used the first estimate of 58m/s as a cutoff, keeping in 
mind that the actual threshold might still vary and is in fact de¬ 
pendent on which of the fragments collide. 

Because of the large number of simulations, we needed an 
automatic filtering algorithm to find the systems that evolved into 
a co-orbital. We used a four-step filtering process. 

1. We selected the systems in which at least two fragments have 
survived until the end of the simulation (20-30 years). Erag- 
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Table 1: Parameters of Janus and Epimetheus (Jacobson et al.||2008 Thomas et al. 2013 i 


parameter 

unit 

Epimetheus 

Janus 

mass 


^epi 

0.11667- 10-*" 

^jan 

0.58333 ■ 10-" 

escape velocity 

m/s 


39 

^esc 

65 

density 

kg/m^ 

Pepi 

640 + 62 

Pjan 

630 + 30 

mean radius 

Ro 


3.85 ■ 10-4 


1 

o 

00 

00 

semi-major axis 

Ro 

^epi 

1.000172 

^jan 

0.999842 

eccentricity 


^epi 

0.0097 

^jan 

0.0068 


Notes, is the mass of Saturn and Rq is the mean distance of Janus and Saturn. 


Table 2: Properties of Saturn as used in the simulations 


zonal gravitational 
coefficients 

J2 

h 

1.62980- 

-9.150- 

10-2 

10-4 



h 

1.030- 

10-4 


mass 


5.688 - 

1026 

kg 

equatorial radius 



0.4 

Ro 

Janus mean distance 

Ro 

1.51464 

- 10*^ 

m 



Fig. 4: Frequency of co-orbitals as a function of the initial speed 

Uinit of the fragments for batch 3 

ments can be lost through ejection, collision with Saturn, or 
collision with each other. 

2. From the remaining systems, we selected the systems with 
exactly one pair of fragments that survived until the end of 
the simulation and overlap in their range of semi-major axes 
for the time after the last collision in the simulation. 

3. From these, we selected the systems in which the mean semi¬ 
major axis of the surviving fragments is smaller than 1.2Ro 
and the minima and maxima of the semimajor axes are no 
more than 0.02R() apart. 

4. We plot the semi-major axes of the remaining fragments 
against time and inspected them visually to check for pe¬ 
riodic crossings. We thus determined false positives of the 
automatic hltering process and removed them. 


3. Simulations and results 

In each simulation, we sampled the parameters Uinit and ip in 
equally spaced steps. We started our investigation with a batch 
of 4000 simulations, covering 20 different initial speeds Uinu and 
200 different angles p (see Tables [^1). The initial speeds were 
between 10 and 100 times as high as the escape velocity of Janus. 
We found that around 1.5% of the simulations evolved into co¬ 
orbital systems. From the analysis of these systems, we were 
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able to establish a common formation scenario as depicted in 
Fig. After the initial collision event, the fragments first spread 
out on orbits around Saturn. Then, they re-collide pairwise, most 
commonly a with y and yS with 5. The velocities of these col¬ 
lision events are very high (lOOOm/s and above), far beyond 
where a perfectly inelastic joining, as is assumed by Mercury, 
is realistic. But the moons created in this way quickly (in just a 
few orbits) enter into co-orbital horseshoe motion after the sec¬ 
ond moon has been formed. No co-orbital systems are the result 
of ejection of fragments or collision with Saturn. 

Then, we investigated whether lower initial velocities (Ta¬ 
bles [^2) result in lower collision velocities. We reduced the pa¬ 
rameter range of Uinit by a full order of magnitude. This indeed 
reduced the collision velocities by nearly an order of magnitude, 
but not to the point where simple merging of bodies was a re¬ 
alistic outcome for any of the collisions. We note, however, that 
co-orbital formation increased to 5.7%. 

Our third batch of simulations was divided into three parts: 
First, we simulated with the same masses as before, but at even 
lower initial velocities (Tables [^3a). As expected, the colli¬ 
sion velocities decrease even more (see Fig.j^for a comparison 
of collision velocities). The portion of co-orbitals increases to 
8.4%. Another result is that the frequency of co-orbital formation 
decreases for very low initial velocities (below ~ 0.015Ro/d or 
around half of the Janus escape velocity). The highest observed 
frequency of co-orbitals is around 0.025R()/d. 

For the second and third part of this batch, we decreased the 
masses of the fragments by a half. One part was ran at initial ve¬ 
locities like in the previous batches (Tables [^3b) and the other 
at lower ones (Tables [^3c). By comparing simulations 3a and 
3c, we see that the mass of the fragments has little effect on the 
distribution of generated co-orbital systems. Only the decrease 
in occurrence of co-orbitals for the lowest initial velocities (see 
Fig.g is less prominent for lower mass. The decrease at these 
initial velocities indicates a lower limit for the formation of co¬ 
orbitals. Supplementary simulations indicate that for even lower 
velocities yinit < 2 ■ 10“^Ro/d, formation of co-orbitals is quite 
unlikely (~ 0.7%) compared to the other simulations. A cut-off 
would probably depend on the mass of the fragments, assum¬ 
ing it is related to the escape velocity of the fragments. This 
could explain the higher formation probability for lower-mass 
fragments at low initial velocities. 

We thus have two constraints on the initial conditions if we 
wish to determine a range likely to spawn co-orbitals: The col¬ 
lision velocities are highly dependent on the initial velocities. 
Fow collision velocities are required for the fragments to col¬ 
lide constructively, thus low initial velocities. At the same time, 
the formation frequency of co-orbitals decreases after a certain 
point, so the initial velocity cannot be arbitrarily low. The exis¬ 
tence and width of a sweet spot for the formation of co-orbitals 
is key to decide the viability of this model. 
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Table 3: Parameters for our simulations 


set 

yinit(Ro/d) 
min max 

steps 

min 

ifi (rad) 
max 

steps 

masses 

integration 
time (yr) 

frequency of 
co-orbitals 

collisions below 
threshold 

1 

0.5 

4.0 

20 

0 

n 

200 

^epi 

^jan 

20 

1.5% 

0 % 

2 

0.05 

0.4 

10 

0 

n 

100 

^epi 

^jan 

20 

5.7% 

0 % 

3a 

0.005 

0.04 

10 

0 

n 

100 

^epi 

^jan 

30 

8.4% 

12.5% 

3b 

0.05 

0.4 

10 

0 

n 

100 

^epi /2 

^^jan/2 

30 

7.4% 

0 % 

3c 

0.005 

0.04 

10 

0 

71 

100 

^^epi/2 

^jan/2 

30 

8.7% 

15.8% 


Notes. For constants, see Table^ 



libration period (J/E libration period) 


Fig. 5; Period of the horseshoe orbits. Note the similar distri¬ 
bution for all masses versus selected masses. The masses of the 
bodies have almost no influence on the period. 


We examined the results of simulation batch 3c more closely. 
In many but not all cases, the fragments collided in such a 
way that the final mass of the bodies was equal to Janus and 
Epimetheus. This is the case for 36 of the 87 co-orbital sys¬ 


tems we obtained in this batch, or 3.6% of all simulations in 
the batch (see Table|^. The collision velocities reach up to about 
450m/s. The collisions occur at angles smaller than 5° for nearly 
all cases. Assuming a limit of 58m/s for constructive collision, 
we And that 15.8% of all collision events in this batch remain 
(see Fig.[^. By multiplying the original success rate of 3.6% by 
15.8% squared (since both collision events need to be below the 
threshold), we obtain a corrected success rate of 0.09% to gen¬ 
erate Janus and Epimetheus or 0.2% to generate any co-orbital. 

We compared the libration period of the horseshoe orbits to 
the period of the horseshoe orbit of Janus and Epimetheus (see 
Fig. 0- Janus and Epimetheus cross orbits roughly every four 
years. This period is highly variable for different simulations 
and shows no particular connection to the mass of the moons. 
Instead, there seems to be a general preference for short periods. 

We see more influence of the mass of the fragments on the 
trajectories by examining the relative range of the semi-major 
axes. With our setup of four fragments, two each of identical 
mass, there are four possible combinations of a Anal co-orbital 
pair, listed in Table 

We now consider the ratio of radial amplitudes 
- max(ai)-min(ai) 

A = -——- . , (3) 

max (a 2 ) - mm (^ 2 ) 
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Table 4: Classes of co-orbital systems by mass. The frequency is relative to all co-orbital systems in simulation batch 3c. Check 
Fig.|^for an example orbit for each class defined here. 


class 

combinations 

mass ratio 

frequency 

A 


{a + y) 


(f3 + 5) 

4.00 

37.5% 

B 

(a) 

iP + y + S) 

(y) 

{a + p + 6) 

0.11 

47.2% 

C 


(a + y + 6) 


{a+p + y) 

1.40 

6.94% 

D 

(a+13) 

{y + 6) 

(a + 6) 

(P + y) 

1.00 

8.3% 



Fig. 6: Four classes of co-orbitals by mass ratio and A (see Ta¬ 
ble |4|lJ[he_variables_m|eirivem proportional for horseshoe or- 
bits ( [Dermott & Murray|198T] Eqn. (1)). The red circle indicates 
the only pair in a tadpole orbit we observed in this batch of sim¬ 
ulations. 

where ai is the set of semi-major axes of the heavier body of the 
hnal co-orbital pair after the last collision event in the simula¬ 
tion, and 02 is the same set for the lighter body. For the Janus- 
Epimetheus system, this is 

Aj/e = 0.2102. (4) 

By plotting A against the ratio of masses (see Eig.[^, we recover 
the expected inversely linear correlation between the two. Each 
class has similar ratio of radial amplitudes. There is one simu¬ 
lation with significant deviance (see Eig. red circle) from this 
law. Closer investigation reveals that in this case, the surviving 
bodies enter into a tadpole orbit (see Eig. [^. We have found no 
other case of tadpole orbits in our simulations. 

4. Conclusion and outlook 

We here investigated the possibility of a formation process of 
the Janus-Epimetheus system starting from a collision event that 
spawns four large fragments. Additional debris was neglected. 
We used numerical integration of the N-body problem to obtain 
the evolution of 8000 initial conhgurations. We selected the con- 
hgurations that led to the formation of a co-orbital system. By 
analysing these, we established a common scenario that leads to 
the formation of co-orbital systems: After the fragments spread 
out on orbit, they re-collide and form a pair of co-orbiting bod¬ 
ies. 
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We investigated the circumstances of these secondary colli¬ 
sions and showed that the initial velocity has a high influence on 
the collision velocities. At the same time, we found that lower 
initial velocities yield a higher probability of co-orbital forma¬ 
tion. This reinforces our assertion that for our assumed initial 
condition to be achieved, the prior collision event must have oc¬ 
curred at a low velocity. It is unclear, however, whether the col¬ 
lision velocities observed in these simulations are low enough 
that a perfectly inelastic collision and joining of the bodies is a 
good approximation. We assume that relative velocities below 
58m/s allow the fragments to collide constructively. Previous 
studies on the collisions of celestial bodies are not fully appli¬ 
cable to our system. Most works dealt either with collisions of 
rocky bodies and/or of a small projectile hitting a large target. 
Therefore it is necessary to study the collision behaviour of large 
low-density fluffy objects such as Janus and Epimetheus in more 
detail. Studies in this area would help to ultimately decide the 
viability of our formation model. We also found that there is a 
lower cut-off below which the formation of co-orbitals decreases 
significantly. We attribute this to the inability of the fragments 
to escape their combined gravity when the kinetic energy is too 
low. This imposes two constraints on the initial velocity: It must 
be low enough to obtain reasonably low collision velocities. At 
the same time, it must be high enough to allow the fragments to 
escape their mutual gravitational pull. Our results indicate that 
this sweet spot is quite narrow, but indeed exists. 

Another limiting feature of our physical model is the number 
of initial fragments adopted in the simulations. To strengthen the 
scenario proposed here, we performed two sets of simulations, 
one considering six, the other eight initial fragments. Consider¬ 
ing initial velocities as those of simulation 3c, 100 simulations 
were performed for each case (six and eight fragments). The re¬ 
sults showed an even higher rate of co-orbital generation (18- 
20%). Therefore, the scenario proposed here is somewhat robust. 

Eurthermore, we are aware that our initial conditions are 
not very natural. Instead, we chose a simple configuration to 
constrain the otherwise quite large number of free parameters. 
The main result of this study is therefore that the proposed 
formation model is possible. We are confident to And co-orbital 
formation for a wider range of initial conditions building 
upon this result. Ideally, a future work would include the initial 
collision as part of the simulation and use a better approximation 
for the secondary collision events. 
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(c) class C (d) class D 

Fig. 7: Example orbits from all four classes. Radial axis is the semi-major axis relative to the mean semi-major axis. Angular is the 
mean longitude in a frame rotating with the mean motion associated with the mean semi-major axis. To view the evolution of each 
of these cases, visit http://goo.gl/JgJ7wZ 
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a + y - /3 + 6 



Fig. 8: Single case of tadpole orbit. The fragments have the mass of Janus and Epimetheus. The plot shows the semi-major axis 
relative to the mean semi-major axis plotted against the mean longitude in a frame rotating with the mean motion of the mean 
semi-major axis. 









